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Abstract 



' The T-noninvariant transverse polarization of neutrons is considered for muon capture by a ® Li 

^ I nucleus with decay into a quasistationary 2"*" state of the three particle a + n + n continuum. This 

• polarization is orthogonal to the plane spanned by the polarization axis of the initial mesic atom and 

' the neutron momentum. The situation in which neutrons that are emitted in the plane orthogonal 

to the axis of the mesic-atom polarization and which have oppositely directed momenta equal in 
magnitude are selected is analyzed. The wave function of the final state is constructed by using the 
method of hyperspherical harmonics. In the approximation of the allowed Gamow- Teller transition 
1+ — > 2^, this neutron polarization is expressed in terms of the T-noninvariant relative phase of 
' reduced matrix elements for transitions from the ^Li ground state to various final-state configurations. 

(N ■ 

■ 1 Introduction 
O' 

' Evidences for violation of time reversal invariance (T-invariance) are presently obtained only in K^-meson 

\ decays. Attempts to find T-noninvariant effects in nuclear reactions, 7-transitions, /3-decays of nuclei and 

free neutrons, along with searches of electric dipole moments of nuclei and elemental particles, have not 
been successful. However, the results of these investigations allowed to lower the upper limit of possible 
mixture of T-noninvariant interactions to ~ 10^"^ — 10^^ (see, e.g., 0). 
^ \ New approaches to T-invariance tests are of obvious interest. In Refs.[^-Q] the possibilities of T- 

C ■ invariance studies in modern experiments on muon capture were considered. In particular, reactions with 

transition to continuum /i + d^n + n + f^ and /i He ^ d + n + were discussed in Ref. . It was 
, - pointed to the fact that neutron polarization along [rifc x n^], where rifc and are unit vectors along 

r> I neutron momentum and muon polarization, respectively, should be sensitive both to T-noninvariant 

■ phases of formfactors in the week semileptonic Hamiltonian and to T-invariance violation in nucleon- 
nucleon interaction. 

Of special interest are few body systems for study of T-noninvariant components of nucleon-nucleon 
potentials. Only in such systems reliable connection between experimental results and constants of 
Hamiltonian may be established. Recently progress was made towards microscopic description oi A — Q 
nuclei as three body systems a -I- N -I- N in the framework of hyperspherical functions method (see, e.g.. 
In particular, in Ref.|^ the calculation was performed of the probability of muon capture by 
nucleus ^Li with transition to the ground state of ^He. 

The binding energy of the nucleus ^He equals 0.975 MeV. Above this value the three particle continuum 
a + n + n begins. There is a resonance in this continuum with spin and parity 2"*", energy 1.8 MeV and 
width 0.1 MeV. Transition into this resonance state as the result of muon capture by nucleus ^Li with 
spin and parity I'*' is of Gamow- Teller type. Thus the existence of the resonance simplify the description 
of muon capture with decay to continuum. This paper is devoted to estimation of T-noninvariant neutron 
polarization along the direction [rifc x n^] after muon capture by ^Li with decay to continuum. 



2 Muon capture — general formalism 

Non-relativistic transition Hamiltonian in nucleon space results from relativistic Hamiltonian by Foldy- 
Wouthuysen transformation (see, e.g., |Q) and is expressed as a power series in the ratio En/Mc^, where 
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is neutrino energy, and M is nucleon mass. To estimate the effect in a Gamow- Teller transition we 
restrict our consideration to zero-order terms. Thus a dimensionless Hamiltonian is of the form 

A 

/i' = ^exp (^-i^^^ (^-«5vb|((T^,cr^) +5iAb+(cr^,(7i,)(Tj-^ -r_(j). (1) 
We use here 4-vcctors 

bA(o'p, CT^) = i Uf,{0, cr^)7A(l + l5)ui,{Pu,cri'), (2) 

which are introduced in accordance with Ref.|^ and expressed in terms of muon and electron 4-spinors. 
The matrixes 7a are taken in pseudoeuclidian metric. The Hamiltonian h' is a matrix in a space of 
spin projections of muon ct^ and neutrino a^. It contains formfactors of vector gy and axial- vector 5^ 
interactions, depending on squared transferred 4- momentum fc^ = k\k\. Summation over all A nucleons 
enters in the Hamiltonian. Spin operator aj and lowering isospin operator T^{j) act in a space of j-th 
nucleon. We have 

f_|p>=|n>, f_|n>=0. (3) 

Nucleon coordinates Vj are reckoned from the center-of-niass of the nucleus. In the non-relativistic 
approximation a muon 4-spinor is of the form 

u.iO,a,)=[^^[^^^), (4) 

where ipfj_{af^) is a usual two-component spinor. 

The nucleus captures a muon from the Is-statc. If nucleus spin Ji differs from zero, the state is splitted 
into two hyperfine sublevels with angular momenta — Ji ±1/2 and energies E^. Interference terms 
~ exp (±i(_E+ — E^)t/fL) of differential capture probability oscillate rapidly and go to zero on averaging 
over life time of mesoatom. So angular correlations in muon capture should be calculated separately for 
each hyperfine sublevel, even though these sublevels are not distinguished in the experiment. 

Let the wave function "^j-Mi describes initial nucleus in a state with projection Mi of spin Ji on an 
axis z. The mesoatom state from which the muon is captured is determined by the function 

|F >=Y.-dF) E (5) 

The amplitudes a^{F) contain an information on mesoatom polarization. They are normalized to the 
unity )P = 1, where ^ is a projection of the angular momentum F on the polarization axis (axis 

z). Let us introduce a unit vector along this axis. Mesoatom polarization is given as usual by the 
quantity 

PliF) = ^^, <e>=5]e|ae(F)p. (6) 

Let the final state of the system a-|-n-|-n, formed after muon capture by the nucleus ^Li, is described 
by the function where cti and 172 are spin projections of 1-st and 2-nd neutrons on an axis z'. 

Thus we have for the amplitude of probability to find the 1-st neutron with spin projection cti and the 
2-nd neutron with spin projection (T2 on an axis z' 

a^(ai,a2) ~ ^a^li^) J] C^^^,,^ < ^f^^\h'\^jai, > . (7) 

The probability of finding of the 2-nd neutron summed over non-observed spin projections of the 1-st 
neutron and neutrino is proportional to the equation 

<-EEi«''(^i'^2)r (8) 

While for the averaged spm projection of the 2-nd neutron we obtain 

<-E^2Ei«''(^i'^2)r (9) 
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Neutron polarization along an axis z' is determined by the ratio 

vl = 4- (10) 

This polarization depends on angular momentum of mesoatom. 



3 Wave functions 

At low excitation energies nuclei with mass number ^ = 6 behave as systems built up from three bodies 
— a-particle and two nucleons. It is convenient to construct the wave functions in the form of series 
on hyperspherical harmonics (see, e.g., [p^). In particular in the framework of this method the wave 
function was calculated of the ^Li ground state 1+ with isospin T = in Ref.@. 

In Ref.[Q the same method was used to study the structure of 2+ continuum state of the system 
a + n + n with isospin T ~ \. Generally a three-body continuum state is described in a center-of-mass 
system by the function ^'p^p,,, depending on asymptotic momenta Pa; and Py, which are conjugated to 
normalized Jacobi coordinates x and y. In the system a + n + n these coordinates 



Y(r„2-r„i), y=y— (r„ ^ I, (11) 

are proportional to the relative radius-vectors of two neutrons p\2 — rn2 ~ r„i and of the a-particle with 
respect to the center-of-mass of the nucleon pair ps = — (r„i -I- r„2)/2. The momenta p^; and Py are 
expressible in terms of neutron and a-particle momenta in the laboratory system 

= ^^(P"2 - Pnl), py = — i= - (p„i +p„2)) • (12) 



/2M^"'" VUd V 2 

Hyperspherical harmonics $'^lm(^) depend on five variables Q — {9, 0XTipx,Oy,ipy) and are fixed by 
five quantum numbers — hypermoment K, orbital momenta l^, ly, total orbital momentum L and its 
projection M on an axis z. Hyperspherical harmonics are of the form 

*^Lm = A^]^'"(sin6')'="(cos6l)'''P^"+^'"+^(cos26') X 



X 



ctZ,^Yi^,nAe.,v.)Yi,,nSey,^y), (13) 



where P"^ are Jacobi polynomials, variable n takes the values n = 0,1,2..., hypermoment K equals 
2n + Ix + ly, and Yim are usual spherical harmonics. Normalized factor is determined by the formula 

M ^ ( 2"!(2» + lx + ly + 2)r(n + -t- + 2) \ ^'^ 

^ I nn + lx + l)V{n + ly + l) J ■ ^ ^ 

Let us introduce wave vectors — Px/h — {kx,Okx,^kx), kj, — Py/h = {ky,6ky,ipky), where polar 
dkx and azimuth Lpi^x angles specify the direction of the vector k^, , as well as the angles 6ky and ipky give 
the direction of ky . Variables k and Ok are determined as follow 



k—^k'^ + k^, kx — ksm0k, ky — kcos9k- (15) 

Note, that the energy in the center-of-mass system equals 

E = (16) 

So parameter 9k determines how this energy is distributed between the subsystems n -I- n and a + (2n). 
In according with Ref. 0| the continuum wave function of a -I- n -I- n may be written as series on the states 
with definite angular momenta J/ and their projections Mf on an axis z 

JfMf K'j nia 
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Here S is the total spin of two neutrons, and the index 7 denotes a set {l^lyLS}. The expansion coefficients 
contain the hyperharmonics, depending on variables flk — {Ok, Okx, Vkx, Oky, ^Pky)- 

Calculations of Ref. have shown that the continuum resonance state 2+ of the a + n + n system is 
formed mainly by two terms of the above expansion, corresponding to the same hypermoment K = 2. 
The contribution of the component 7I = {0220} in the internal region varies from 45 to 70 % depending 
on the calculation method, while the contribution of the component 72 = {1111} changes from 35 to 
20 %. We shall restrict our consideration to these two terms in the expansion (|l^) for the wave function 
of the final state. 

The wave functions ^'j^M; and ^'pjp^ are, of course, antisymmetrized over transposition of space, 
spin and isospin coordinates of two nucleons of the system a + N + N. Therefore matrix element of 
Hamiltonian h' = h'l + h'2 is equal to the sum of two identical matrix elements of the operators h'l and 
/ij. In this paper we take into account only the contribution of the allowed s-wave matrix element to the 
amplitude of Gamow- Teller transition 1+ — s- 2^. Introducing the reduced matrix elements in accordance 
with Ref . Ipl]] , we get for the matrix element of spherical component of the operator a 

< *J;m>oI*^.m. >^ ^c'/jj%m]K,. (18) 

Spherical components are determined by usual rule a±i = =f('5'x i ay)/\/2), ao = Oz- 

If T-invariance holds, thus the phases of the wave functions may be chosen in accordance with the 
standard condition 

flJM >= (-1)^+*^|J-M >, (19) 

where T is the time reversal operator. In this case the reduced matrix elements [lOlJif-^ are real. In 
another way the coefficients of expansion of the wave functions j^Mi and "^^^Mf '-'^ hyperspherical 
harmonics are hyperradial functions which satisfy the set of coupled differential equations of second 
order. The reduced matrix elements are the one-dimension integrals of these functions. Following the 
rule (|l9|), one may make real all hyperradial functions and, consequently, reduced matrix elements. 

With nucleon-nucleon potentials violating time reversal invariance the set of coupled differential equa- 
tions have no real solutions. The rule (p^ ) does not hold. So the reduced matrix elements become complex. 
The procedure of calculating of imaginary T-noninvariant corrections to functions satisfying the set of 
coupled differential equations of second order was considered in Ref. [|l3j . In this paper we express the 
neutron polarization induced by time reversal violation in terms of imaginary parts of matrix elements. 



4 Analysis of T-noninvariant effect 

In the framework of our formalism the neutron polarization (|l^) depends on three momenta — on neutrino 
momentum p^, which is equal by magnitude and opposite by direction to the momentum of the center- 
of-mass of the system a + n + n, and on Jacobi momenta p^; and Pj,. In principle, these three momenta 
may be determined by measuring the momenta Pni, p„2 and in the laboratory coordinate system. 
However, such experiment is non-realistic. 

As the other limiting case one may consider the situation when the momentum of only one neutron is 
measured, while the integration over all possible momenta of the other neutron, a-particle and neutrino 
holds. In such experiment it is impossible to control the decay into the resonance state 2+. Therefore 
this case is not interesting for us. Note in addition, that recalculation of the effect, expressed in terms 
of natural for three-body problem Jacobi momenta p^,, p^, to laboratory momenta Pni, p„2 and p^ 
needs the cumbersome numerical procedure. Such recalculation would complicate the interpretation of 
observable effect. 

The optimal experiment should be not-too-compex on setting, on the one hand, and simple on in- 
terpretation, on the other hand. We consider the following situation. Lets assume that the momenta 
of two neutrons are detected in the experiment. We select the events when these momenta are equal 
by magnitude and opposite by direction. Thus the magnitude and direction of Jacobi momentum p^, is 
determined by equation (l2h. Furthermore, the total momentum of two neutrons equals zero, therefore in 
accordance with formula (|12|) the momentum pj^ coincides by direction with the a-particle momentum. 
Besides, due to the fact that the total momentum of four particles in the final state equals zero 

Pnl + P«2 + Pq + Pi. = 0, (20) 
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in the selecting cases Pa = — Piy- This means that the direction of Jacobi momentum Pj, is exactly 
opposite to the neutrino momentum. Thus the integration over aU non-observable direction Py and Pi, 
easily performs. To control the decay into the resonance 2"*" continuum state one need to measure only the 
total energy of the a-particle. In addition the integration over the parameter 9k holds if the distribution 
of the total energy between the subsystems n + n and a + (2n) is not taken into account. 

We normalize the continuum wave functions of the system a + n + n as in Ref . . Thus we obtain 
the following expression for the differential probability (||) of detecting of two neutrons with opposite 
momenta (Pni = —Pn2) and with Jacobi momentum p^; in the solid angle dflx for the Gamow- Teller 
transition Ji ^ J f = 2 

dwE = A,\gAf,^j±\c{J.,F) (|[101],ip + |[101],.p) ^^rffc. (21) 
Here we use the constant 

^ 3 il + EjEf){l + m^cyE,r' ^^^> 

where 



G is the week constant, 0c is the Cabbibo angle, is the muon mass, Z is the charge of the initial 
nucleus, and R{Z) is the correction for nucleus size Neutrino energy E^, is determined by the total 
energy Ef of the final system a + n + n, which includes rest masses of the particles, and by the total 
energy Q^, released in muon capture, as follow 



E^ = Ef 



2Q^ 
Ef 



1/2 



Q,[l-^^+...]. (24) 



Differential dk is related with differential dE of the final system energy according the formula ([iq): 
dk = dE/{h^/2E), where E and dE should be taken equal to the energy and width of the resonance 
state 2+. The solid angle Afiy fixes the value of the possible deviation of the momentum pj, from the 
direction, opposite to neutrino momentum. The coefficient C{Ji, F) depends on spin Ji of initial nucleus 
and on total angular momentum of mesoatom = J; ± 1/2 and is of the form 

C( J„ F) = 1 + V6U{F^Jd, J.i)C/(2Jai, IJ,), (25) 

where U{abcd, ef) = ((2e + 1)(2/ + l)y/'^W(abcd, ef) is the normalized Racah function @. 

Note, that the differential probability ( pl|) does not contain the terms like ~ cos 6*, where 9 is the 
angle between the momentum p^; and the mesoatom polarization axis n^. This is associated with the 
identity of emitting neutrons and, consequently, with ambiguity of the choice of the direction p^;. There 
exists also the other special feature of the coefficient C{Ji = 1,F). Let Ji — 1, thus F = 1/2 or 3/2. 
It is easy to check that C{Ji — 1,F) equals zero, if the total angular momentum of the initial state is 



F = 1/2. This is the result of angular momentum conservation. Indeed, the expression (21) is obtained 
in the approximation of the s-wave neutrino emission with respect to the center-of-mass of the system 
a + n + n. It is clear that as neutrino spin equals 1/2 and the spin of the final resonance state is Jf — 2, 
the total angular momentum in the final state is 3/2 or 5/2. So this state can not be obtained from the 
initial state with F — 1/2. 

We turn now to the situation when neutron momenta p„i = — p„2 are perpendicular to the polar- 
ization axis of the initial mesoatom. By convention we assign a number 1 to one neutron and a 
number 2 to the other neutron and fix the direction of the vector p^;. According to the equation ( p^ this 
vector coincides by direction with momentum p„2- Choosing an axis z' along [p^; x n^], we calculate the 
polarization of the 2-nd neutron along the axis z' using the formulas (|)-(|l^). We get 

^_ 32 [6 DiJ.,F) Im ([101],i[101];,) 
^1 2771 V 5^'^ ^C(J„f^) |[101]^i|2 -I- |[101]^2P- ^ ^ 
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Under a different choice of neutron numbering we change only the direction of the vector p^; and, con- 
sequently, of the axis z' . Thus we should interpret the result ( p6| ) as polarization of any of two identical 
neutrons along [n/j x n^], where is the unit vector along the neutron momentum. The expression 
obtained contains the spin factor 



D{J,,F) 



i2J, + 1)F 
F+1 



1/2 



{2E+1)U{121E,12) X 



E=0,2 



Ui-J,F1,FJ, 




J^ F 1/2 

J, F 1/2 )U{lJaE,2J{) 
E 1 1 



(2f +1)(2J, + 1) 
3 



1/2 



E (2G+1) 



1/2 



G=0,2 



H=0,2 



'2H +l)U{nHE,Gl) 



F 


1/2 




F 


1/2 


I- 




1 






J^ 








ill 




G 





(27) 



It equals zero when Jj — 1, F — 1/2 for the same reason as C{Ji,F). For = 1, F = 3/2 we have 
C(l,3/2) = 1.5, 1^(1,3/2) = 0.237. 

We see that the neutron polarization ( p6| ) differs from zero if the phases of the reduced matrix elements 
[lOlJ^i and [101]-y2 do not coincide. Let us remind that if T-invariance holds the phases of wave functions 
may be chosen according the rule (|19|). In this case all reduced matrix elements are real. Denoting 



T-noninvariant phase of the product [101]-yi[101]*2 by 0, where according Ref. (f) < 10 we get a 
rough estimate of the effect 



O.10Pi(F) 



|[101],2/[101]7ll 
l + |[101]^2/[101]^l| 



(28) 



Note, that the polarization pi{F) of mesoatom /i -|-^ Li in the state F — 3/2 is about 0.05 ||T^. How- 
ever, in the experiments with previously polarized atoms (see, e.g., jl^, 0) the quantity Pi{F) may be 
considerably higher. 



5 Conclusion 

In this paper T-noninvariant transverse polarization is considered for neutrons from reaction of muon 
capture by ^Li nucleus with decay to resonance 2+ continuum state of three particles a -I- n -I- n. This 
polarization is normal to the plane formed by polarization axis of initial mesoatom and neutron momen- 
tum. We analyze the situation when selected neutrons have equal and opposite directed momenta which 
are perpendicular to the axis of mesoatom polarization. 

To estimate T-noninvariant effect we use an explicit expression for the wave function of the final 
state in the hyperspherical harmonics method. An analysis of the structure of 2+ resonance, performed 
in Ref.j7|, allowed us to take into account only two leading terms in this wave function. Only the 
allowed matrix element of Gamow- Teller transition 1+ 2+ was considered. In these approximations 
the neutron polarization is expressed in terms of T-noninvariant relative phase of reduced matrix elements 
for transitions from the ground state of ^Li nucleus into the different configurations of the final state. 

The current level of microscopic description of nuclei A = 6 as three-body systems a -|- N + N in the 
method of hyperspherical harmonics allow to calculate the reduced matrix elements. We intend to 

perform such calculations as the next step. This will enable to relate the effect considered with parameters 
of nucleon-nucleon potentials violating time reversal invariance (see, e.g., |]l8|). 

The author is grateful to Yu.V.Gaponov, B.V.Danilin and N.B.Shul'gina for useful discussions. This 
work was supported by International Science Foundation, grant number is M7C300. 
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